Anisotropic singularities in modified gravity models 
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We show that the common singularities present in generic modified gravity models governed by 



actions of the type S = Jd Xy/—gf(R, 4>, X), with X 



jg a d a 4>db4>, are essentially the same 



anisotropic instabilities associated to the hypersurface F((j>) — in the case of a non-minimal 
coupling of the type F{cj>)R, enlightening the physical origin of such singularities that typically arise 
in rather complex and cumbersome inhomogeneous perturbation analyses. We show, moreover, that 
such anisotropic instabilities typically give rise to dynamically unavoidable singularities, precluding 
completely the possibility of having physically viable models for which the hypersurface ^ = is 
attained. Some examples are explicitly discussed. 



PACS numbers: 98.80.Cq,98.80. Jk,95.36.- 



I. INTRODUCTION 

In the absence of a more fundamental physical model 
based on first-principles for the description of the cos- 
mic acceleration discovered more than a decade ago[lJ 
(see, for reviews, [2|), many dark energy phenomcno- 
logical models have been proposed and investigated in 
detail. In particular, the questions about the stability 
against small perturbations in the initial conditions and 
in the model parameters are always the first requirement 
demanded to assure the physical viability of any cosmo- 
logical model. The most part of such dark energy models 
belong to the general class of cosmological models gov- 
erned by an action of the type (see, for instance, [3]) 



S = 



(1) 



where R stands for the spacetime scalar curvature, is a 
scalar field, X = — T;g ab d a d)db4>, and / is a smooth func- 
tion. Quintessence models|J], for instance, correspond to 
the choice f(R, <j), X) = j±-R-±g ab d a 4>d b <j) + V((t)). Non- 
minimally coupled models [a, 0,0, on the other hand, are 
typically of the type 

f(R, 0, X) = F{4>)R - ^g ab d a <t>d b <l> + V{ct>). (2) 

Many other models discussed in the literature correspond 
yet to the case f(R,(f>,X) — g(R,<f>) + h((j),X), includ- 
ing fc-essence|8] and the string- inspired case of a Dirac- 
Born-Infeld tachyonic action[9(. (For more recent works, 
see [13].) The particular case of pure modified grav- 
ity f(R,4>,X) = f(R) (see, for a recent review, [ill ]) 
has been intensively investigated as an alternative to 
quintessence. Some primordial inflationary models [l2T] 



are also described by actions of the type ([I]) . Since one of 
the proposals of any cosmological model is to describe our 
universe without finely-tuned parameters, a given dark 
energy or inflationary model would be physically viable 
only if it is robust against small perturbations in the ini- 
tial conditions and in the model parameters. This is the 
question to be addressed here. 

Non-minimally coupled models of the type ([2]) are 
known to be plagued by anisotropic singularities in the 
phase space region corresponding to F((f>) — 0. For in- 
stance, Starobinski[l3j was the first to identify the sin- 
gularity corresponding to the hypersurfaces F{4>) = 0, 
for the case of conformally coupled anisotropic solu- 
tions. Futamase and co-workers[5j identified the same 
kind of singularity in the context of chaotic inflation in 
F{4>) = l-£</> 2 theories (See also [3]). In [3, it is shown 
that such kind of singularities are generically related to 
anisotropic instabilities. 

Many authors have described different singularities 



corresponding to ^ — in general models like (|TJ) (see, 
for instance, [3() or, more commonly, in pure f(R) gravity 
models (see, for instance, [Hl[l3])- Such singularities ap- 
pear typically in rather complex and cumbersome inho- 
mogeneous perturbation analyses, obscuring their phys- 
ical origin and cause. In this work, we show that these 
singularities are essentially due to anisotropic instabili- 
ties, in a similar way to those ones described in [l4j for 
models of the type @. Moreover, we show that such in- 
stabilities typically give rise to dynamically unavoidable 
singularities, rendering the original model physically un- 
viable. 

One can advance that there are some geometrically 
special regions on the phase space of the model in ques- 
tion by an elementary analysis of the equations derived 
from the action ((T|). They are the generalized Klein- 
Gordon equation 
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D a (f,xd a (f>) + U = 0, 



(3) 



2 



and the Einstein equations 



1 



FG a b = -{f-RF)g ab + D a D b F-g ab DF 
- -f,xd a <t>db<i>, 



(4) 



where F — F(R,<p,X) = We will consider here the 
simplest anisotropic homogeneous cosmological model, 
the Bianchi type I, whose spatially flat metric is given 
by 

ds 2 = -dt 2 + a 2 {t)dx 2 + b 2 {t)dy 2 + c 2 (t)dz 2 . (5) 
The dynamically relevant quantities in this case are 

b 



a b c 

Hi = -, H 2 = t, and H 3 = - 
a b c 



(6) 



For such a metric and with a homogeneous scalar field 
4> = </>(£), Einstein Eq. (H} can be written as 



(7) 



FG m = ~{f-FR)-(H x +H 2 +H 3 )F 
FC 1 

—T" = ^(f-FR) + (H 2 +H 3 )F + F, (8) 
FC 1 

IJJ = -(f-FR) + (H 1+ H 3 )F + F, (9) 



b 2 
FG 33 



(/ - FR) + (H 1 + H 2 )P + P, (10) 



c 2 2 

and the generalized Klein-Gordon equation will read 
j t (f, x fy + {H x +H 2 + H 3 )f tX j> - /, = 0. (11) 

Notice that (fTTj) is a second order differential equation 
for (j>, while Eqs. (fT0|) form a higher order system of 
ordinary differential equations. Since F — F(R,<fi,X), 
the term corresponding to F involves, in fact, sec- 
ond derivatives of R and, consequently, third deriva- 
tives of Hi, i = 1,2,3. Thus, the corresponding phase 
space M. is 11-dimensional and spanned by the variables 
{^H X ,H X ,H X ,H 2 ,H 2 ,H 2 ,H 3 ,H 3 ,H 3 ). Eq. © corre- 
sponds to the energy constraint. It restricts the solutions 
of (|5))- (fTTj) on a certain (vanishing energy) hypersurface 
£ of M. Thus, effectively, the solutions of ([^- (TTTj) are 
constrained to the 10-dimensional manifold £ e M. 

It is quite simple to show that Eqs. (|5j)- (fTU|) are not 
compatible, in general, on the hypersurface T of M cor- 
responding to the region where F(R, <fi, X) = 0. Sub- 
tracting © and (HHJ) from © we have, respectively, on 
such hypersurface 

(H x - H 2 )F = 0, and (H x - H 3 )F = 0. (12) 

Hence, Eqs. ([5|)- (ITU|) cannot be fulfilled in general for 
anisotropic metrics. As it will be shown, the hypersur- 
face J- indeed corresponds a geometrical singularity for 



anisotropic spacetimes which cannot be dynamically pre- 
vented in general by requiring, for instance, that F = 
on the hypersurface F as suggested naively from (TT2"j) . 
Furthermore, the Cauchy problem for the Eqs. (|5j)- (fTTj) 
is ill-posed on this hypersurface, since one cannot choose 
general initial conditions on it. 



II. THE SINGULARITY 

In order to study the geometrical nature of the singular 
hypersurface T, let us consider the Einstein Eqs. ([7])- 
dTTJJ) in detail. For the metric we have the following 
identities 



Goo = H X H 2 + H 2 H 3 + H X H 3 , 



(13) 



Gn = a 1 [H x + H X (H X +H 2 + H 3 )--R) , (14) 



G 22 = b z [H 2 + H 2 (H X + H 2 + H 3 ) - -R ) , (15) 
G33 = c z ( H 3 + H 3 {H X +H 2 + H 3 ) -\r), (16) 



(17) 



R = 2{H X +H 2 +H 3 + H 2 +H 2 + H 2 
+H X H 2 + H 2 H 3 + H X H 3 ). 



Now, we introduce the new dynamical variables p = H x - 
H 2 + H 3 , q = H x - H 2 , and r = H x - H 3 . Notice that 



R = 2p+- (2p 2 + q 2 +r 2 - qr) 
o 



(18) 



implying that R involves terms up to third order deriva- 
tive in p and up to second order in q and r. In terms of 
the new dynamical variables, Einstein Eqs. ([8j)- (fT0|) can 
be cast in the form 

3F = (p + p 2 )F-^f-2 P F, (19) 

qF = -(q + qp)F, (20) 
rF = -{r + rp)F. (21) 

As to the energy constraint we have 

- (p 2 + qr -q 2 -r 2 )F + pP 

+ \(f-FR) = ^l x 4> 2 , (22) 

and the generalized Klein-Gordon equation (J3]) reads sim- 
ply 

(f,X + /WO") h(f,XRR + f,X*<j> + I' f x) 4-U = 

(23) 

Notice that the Eqs. (fT9]) - ([23)) do not involve the terms 
q and r. Moreover, Eqs. (|2ll|) and (|2"l"T) are, respec- 
tively, first order differential equations for q and r, from 
which the terms involving first and second derivative of 



3 



q and r present in the terms F and F of (|T9|) and (|22[) 
can be evaluated directly. The order reduction of the 
system of differential equations attained with the intro- 
duction of the new dynamical variables implies that the 
phase space M. is not 1 1 , but 7-dimcnsional and spanned 
by the variables (</>, 4>,p,p, p, q, r). The solutions are still 
constrained to the hypersurface £ € M. corresponding to 
the energy constraint (|2"2")l . It is clear, however, that the 
manifold £ is, in fact, 6-dimensional. 

There is still a further dynamical restriction on the 
solutions of (QUI-©. From ^ and fH), one has 



rq — qr = 0, 



(24) 



implying that q(t)/r(t) is a constant of motion fixed only 
by the initial conditions. Suppose the initial ratio is 
q(0) / r (0) = 7: this would imply that (H 1 -H 2 ) = "/{Hi - 
H 3 ) for all t, leading to, for instance, c 7 (i) cx a 7-1 (£)&(£) 
in the metric ([5]) . This simplification is a consequence of 
the scalar character of our homogeneous source field, and 
it is also present LJ| in the non-minimally coupled case 
given by actions of the form @ . Let Q be the hypersur- 
face corresponding to q/r constant. Finally, the solutions 
of (| 1 9 [) - ([23]) are necessary restricted to the 5-dimensional 
submanifold Q fl £ of M.. 

A closer analysis of Eqs. ([2H)l and (f2"Tj) reveals the 
presence of the singularity. They can be written as 



* = - ( p+ f) T - 



(25) 
(26) 



In general, the right-hand side of these equations diverge 
on the hypersurface F corresponding to F{R, <f), X) = 0, 
unless q = r = 0. The first observation is that ([25]) 
and (|26p imply, in general, that, if q (or r) vanishes for 
some t, it will vanish for any t. This is why such kind of 
singularity can be evaded in homogeneous and isotropic 
situations. We will return to this point in the next sec- 
tion, with an explicit example. For any physically viable 
cosmological model, small amounts of anisotropy, corre- 
sponding to small q and r, must stay bounded during 
the cosmological history. In fact, it is desirable that they 
diminish, tending towards an isotropic situation. How- 
ever, this does not happen in general if F(R, <fi, X) = 
in ([23]) and (|2l)|) . Let us assume that F ^ on the hy- 
persurface F. (We will return to this point latter.) In 
this case, if any anisotropic solution crosses F , necessar- 
ily q and r will diverge, corresponding to a real spacetime 
geometrical singularity, as one can check by considering 
the Kretschman invariant I = R a bcd,R abcd , which for the 
metric ([5]) is given by 



i 7 



Hi + H{ 



H\Hl 



H l H 3 



H 2 + H; 

tj2 tt2 



H 3 + Hi 



As one can see, I is the sum of non negative terms. More- 
over, any divergence of the variables Hi, H 2 , H 3 , or of 
their time derivatives, would suppose a divergence in /, 
characterizing a real geometrical singularity. Since the 
relation between the variables p, q, r, and Hi, H 2 , H 3 is 
linear, any divergence of the first, or of their time deriva- 
tive, will suppose a divergence in I. 

There are two basically distinct situations where the 
singularity corresponding to the hypersurface T could 
be evaded dynamically. We will show that both are very 
unlike to occur in physical situations. The first one cor- 
responds to the case when the hypersurface T belongs to 
some dynamically unaccessible region. In such a case we, 
of course, do not face any singularity, since F(R,<p,X) 
will never vanish along a solution of the system. This 
would be equivalent to state that T n £ fl Q = 0. In our 
case, it would imply, from (|22|) . that the equation 



(28) 



P F - Xf. x - l -f 



has no solution in M. This would correspond to a quite 
concocted and artificial function /. In particular, for all 
models we could find in the literature having the hyper- 
surface T , the equation (|2"8")) has solutions. 

The second situation corresponds to the already men- 
tioned case where F = on the hypersurface T . From 
(|28[). we see that this requires necessarily that the func- 
tion / be homogeneous of degree ^ in the variable X on 
T . Again, a highly artificial situation. 

Any point on the energy constraint hypersurface £ is, 
in principle, a dynamically possible point. Moreover, it is 
desirable for any cosmological model free of finely-tuned 
parameters that any point or, at least, a large region of 
£ could be chosen as the initial condition for a cosmolog- 
ical evolution. This, of course, includes also the neigh- 
borhood of the hypersurface F provided that F D £ ^ 0. 



III. AN EXPLICIT EXAMPLE 

The singularities described in the precedent section im- 
ply that any model governed by an action of the type 
([1]) having a hypersurface F will certainly present severe 
anisotropic instabilities that will render it physically un- 
viable. Let us work out an explicit example in order to 
illustrate the dynamical role of such anisotropic instabil- 
ities. The pure modified gravity model 



f(R) =R-aR*1n.[l 



R 

R* 



(29) 



where a and i?» are free positive parameters, was recently 
proposed[l6| as a viable model to describe the recent cos- 
mic acceleration. Such a model has a hypersurface F 
corresponding to f'(R) = F(R) — 0, where 



(27) 



F(R) = 1 



(30) 
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In [16| , it is assumed a universe rilled with radiation and 
dark matter, but, for our purposes here, it is enough to 
consider the pure geometrical Lagrangian given by ([29]) . 

Let us start, as in [lfjj], assuming a homogeneous and 
isotropic universe Hi — if 2 = if 3 = if. Einstein Eqs. 
([7|)- (fT0|) for this case would correspond simply to the en- 
ergy constraint 

6HRF'(R) + f(R) - RF(R) + 6H 2 F(R) = 0, (31) 

and to the generalized Friedman equation 

RF'(R) + (2HF'(R) + RF"(R)^ R (32) 

+ lf(R)-(H + 3H 2 ^F(R) = 0, 

where R = 6 if +12H 2 in this homogeneous and isotropic 
case. Note that 

F'(R) = aR * „ > 0, (33) 
(R + R*f 

for R + R* ^ 0. Eq. |j32|) is a third order differential 
equation for if. Hence, the relevant phase space is 3- 
dimensional and spanned by the variables (if, if, if), but 
the solutions are in fact constrained to the 2-dimensional 
manifold £ corresponding to the energy constraint (|31[) . 
The manifold £ is an ordinary smooth surface, with a 
single value of if assigned to each pair (if, if) , provided 
if ^ and R + R* ^ 0. Thus, the solutions of (J32J) can 
be conveniently projected on the plane (H,H), without 
any loss of dynamical information. 

It is convenient to work with the dimensionless quan- 
tities if = ^/R^h, y/R^T = t,R = pi?*. The phase space 
for this model is quite simple. There are only two fixed 
points corresponding to h = ±y/p/12, where p is the 
positive solution of the equation 

2aln(l+p) =P + "Y+ — ' ( 34 ) 

This solution exists and is unique provided that a > 1. 
Both equations (|3"Tj) and (f3"2"| are invariant under the 
transformation r — > —t and h — ► —h, implying that the 
h negative portion of the phase space can be obtained 
from the positive one by means of a time reversal opera- 
tion. Typical trajectories projected on the (h, h) plane of 
the phase space are depicted in Fig. [1] Notice that the 
surfaces corresponding to p = R/R* constant are simple 
parabolas 6h + I2h 2 = p in the (h, h) plane. For the 
model in question, the T surface corresponds to one of 
these parabolas, namely p — a — 1. The point here is 
that the existence of such a surface does not imply any 
singular behavior for the equation l|32p . For instance, 
the solutions crossing the surface F(R) — depicted in[T] 
are perfectly regular. Hence, homogeneous and isotropic 
solutions can cross without problem the singular hyper- 
surface. 

Suppose now the system has a small amount of 
anisotropy, i.e., \q\ <§; \p\ and \r\ <§; \p\. In this case, 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
h 



FIG. 1: Typical trajectories for the model (|29[) projected on 
the (h, h) plane, with a = 2 (the same value adopted in the 
analysis of [IB]). The (attractive de Sitter) fixed point in this 
case corresponds to h w 0.7157. The p-constant parabolas are 
shown. The solid one, in particular, correspond to p — a — 1 
(the singular T surface, defined by F(R) = 0). The regions 
corresponding to F(R) > and F(R) < are, respectively, 
the region above and the region below the parabola p = a— 1. 
As one can see, homogeneous and isotropic solutions can cross 
without problems the singular surface. 

we have from (fl8|) R « 2p + |p 2 and the equation (fl~9|) 
for p will be essentially the same (|3"2"|) obeyed by if in 
the isotropic case, provided the anisotropy is indeed kept 
small along the solutions. For the amounts of anisotropy 
q and r, however, the relevant equations will be (|25[) and 
(l2l))) . Since we know from (|2"4"|) that r(t) — jq(t), we can 
consider here only the variable q 

i= -(* + ?ci A )* (35) 

In any region of the phase space far from the surface 
F(R) = 0, the right-handed side of ()35[) is well behaved. 
Moreover, from the energy constraint (I22|) . we have 

F'(R) ■ 2 7 2 - 7 +l g 2 1 / f(R)\ 

V H ; — rR = — V H -K ; — r 

y F(R) X 3 p 2p \ F(R) J 

(36) 

A closer analysis reveals that 

RF(R) - f(R) = aR, (hi(l + p) - ^-J > 0, (37) 

with the equality holding only for p = 0, implying that 
for the region F(R) > 0, at least, the quantity be- 
tween parenthesis in (|35p is positive, leading indeed to 
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an isotropization of the solutions. For regions close to 
the surface F(R) = 0, on the other hand, the situation is 
qualitatively different. From (|36|) . we have that the right- 
handed side of (|35[) diverges on the surface F(R) = 0. 
If an anisotropic solution reaches such surface, we have 
from Eq. (|35|) that q diverges, implying that this model 
does not admit any amount of anisotropy at all, preclud- 
ing any possibility of constructing a realistic model based 
solely in the geometric Lagrangian (f29|) . Similar results 
hold also for the other functions f(R) discussed in [l6| , 
namely 

f(R) = R-aRjl + — \ , (38) 
with (3 £ (0,1). 

IV. FINAL REMARKS 

The singularities associated with the hypersurface 
F(R,(j),X) = described here are not new. They have 
been discovered and rediscovered many times for many 
different models in rather complex and cumbersome in- 
homogeneous perturbation analysis around a given well 



behaved background solution. Our results, however, en- 
lighten the physical origin of such singularities. They 
arise already in the background level and are related to 
anisotropic expansion rates. Any solution crossing the 
hypersurface F(R,<f>, X) — will not admit, in general, 
any amount of anisotropy, otherwise it will certainly de- 
velop a catastrophic geometrical singularity with, for in- 
stance, the blowing up of the Kretschman invariant (|27[) . 
This, in fact, precludes the possibility of constructing a 
realistic model with solutions crossing the hypersurface 
F(R, <j>, X) = since we would have qualitatively distinct 
behavior for arbitrarily close homogeneous solutions: a 
perfect isotropic and a slightly anisotropic one. 
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